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ABSTRACT. We apply a novel one-dimensional glacier hydrology model that calculates hydraulic head 
to the tidewater-terminating Sermeq Avannarleq flowline of the Greenland ice sheet. Within a plausible 
parameter space, the model achieves a quasi-steady-state annual cycle in which hydraulic head 
oscillates close to flotation throughout the ablation zone. Flotation is briefly achieved during the 
summer melt season along a 7 km stretch of the ~50 km of flowline within the ablation zone. Beneath 

the majority of the flowline, subglacial conduit storage 'closes' (i.e. obtains minimum radius) during the 
winter and 'opens' (i.e. obtains maximum radius) during the summer. Along certain stretches of the 
flowline, the model predicts that subglacial conduit storage remains open throughout the year. A 
calculated mean glacier water residence time of ~2.2 years implies that significant amounts of water are 
stored in the glacier throughout the year. We interpret this residence time as being indicative of the 
timescale over which the glacier hydrologic system is capable of adjusting to external surface meltwater 
forcings. Based on in situ ice velocity observations, we suggest that the summer speed-up event 
generally corresponds to conditions of increasing hydraulic head during inefficient subglacial drainage. 
Conversely, the slowdown during fall generally corresponds to conditions of decreasing hydraulic head 
during efficient subglacial drainage. 


1. INTRODUCTION 

In situ and remotely sensed observations (Zwally and others, 
2002; Joughin and others, 2008; Shepherd and others, 2009; 
Bartholomew and others, 2010) have demonstrated that the 
surface velocity of marginal ice in West Greenland exhibits 
an annual velocity cycle, with peak velocities occurring 
during the summer melt season. Increased summer ice 
velocities have been interpreted to reflect enhanced basal 
sliding due to increased delivery of surface meltwater to the 
bed. Episodic supraglacial lake drainage events, which have 
been inferred to temporarily increase subglacial water 
pressure and result in brief high-velocity events, may be 
overlaid on this annual ice velocity cycle (Box and Ski, 
2007; Das and others, 2008). While the displacement 
associated with surface meltwater-induced acceleration is 
a significant fraction of annual displacement in relatively 
small land-terminating glaciers such as Russell Glacier near 
Kangerlussuaq, it is a much smaller fraction of annual 
displacement in major marine-terminating outlet glaciers, 
such as Jakobshavn Isbrae near llulissat (Joughin and others, 
2008). In marine-terminating glaciers, the primary seasonal 
driver of ice velocity is believed to be the annual advance 
and retreat of the tidewater terminus, which modulates 


glacier flow through a seasonal back-stress cycle (Meier and 
others, 1994; Vieli and others, 2000; Howat and others, 
2005; Joughin and others, 2008). Projections of increased 
Greenland ice sheet surface meltwater production over the 
next century (Hanna and others, 2005) provide an impetus 
to understand the physical basis of the annual velocity cycle 
of the Greenland ice sheet. Recent theoretical (Schoof, 
2011) and observational (Sundal and others, 2011) studies 
of the Greenland ice sheet suggest that future increases in 
surface meltwater production may result in a transition to 
more efficient subglacial drainage and a net decrease in 
basal sliding velocity. 

Studies of alpine glaciers suggest that changes in basal 
sliding velocity are due to a combination of changes in the 
rate of glacier water storage (i.e. total glacier water input 
minus output; Fountain and Walder, 1998; Anderson and 
others, 2004; Bartholomaus and others, 2008) and changes 
in flotation fraction (the ratio of subglacial water pressure to 
basal ice pressure; Iken and others, 1983; Kamb and others, 
1 994). Thus, both the total amount of water storage, which is 
related to englacial water-table elevation, and its rate of 
change influence basal sliding velocity. This explains why 
'bursts' of basal motion are associated with meltwater 
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Fig. 1. The Sermeq Avannarleq flowline (black line) overlaid on 
winter 2005/06 interferometric synthetic aperture radar (InSAR) 
ice surface velocities (Joughin and others, 2010). The locations of 
JAR2 station and Swiss Camp are denoted with stars. The dashed 
line represents the shortest path to the ice-sheet margin from 
Swiss Camp. 


'pulses', while sustained meltwater input, which eventually 
leads to the establishment of efficient subglacial conduits 
and a negative rate of change of glacier water storage, does 
not lead to sustained basal sliding. Changes in the rate of 
glacier water storage, dS/df, are due to changes in both the 
rate of meltwater production (i.e. glacier 'input') and the rate 
of water loss from a glacier, governed by the efficiency of 
subglacial transmissivity via cavities and conduits (i.e. 
glacier 'output'). At the onset of the melt season, the initial 
surface meltwater input exceeds the transmissivity of the 
nascent subglacial system. This increases pressure in sub- 
glacial cavities and conduits, enhancing basal sliding (e.g. 
Fountain and Walder, 1998; Anderson and others, 2004; 
Bartholomaus and others, 2008). Following the onset of 
melt, both meltwater input and subglacial transmissivity 
generally increase with time, the former a direct response to 
meteorological forcing and the latter a response to the 
widening of the ice-walled conduits due to the dissipation of 
energy created by viscous friction (Rothlisberger, 1 972; Nye, 
1976). Enhanced basal sliding is maintained as long as 
meltwater input exceeds subglacial transmissivity (i.e. dS/ 
df>0, or increasing hydraulic head) and is terminated once 
subglacial transmissivity exceeds dwindling meltwater input 
(i.e. d5/df<0, or decreasing hydraulic head). A glacier may 
experience internal meltwater generation due to geother- 
mal, deformational and frictional heating throughout the 
year. While subglacial cavities and conduits typically close 
or diminish in size through internal deformation following 
the melt season (Hock and Hooke, 1993), basal sliding may 
occur throughout the winter as long as internal meltwater 
generation exceeds subglacial transmissivity. 

The desire to quantify the hydrologic contribution to 
enhanced summer basal sliding in the Greenland ice sheet 
provides an impetus to model the magnitude and temporal 
distribution of changes in d5/df within the ice sheet. In this 
paper, we develop a one-dimensional (1-D) hydrology 
model to investigate the annual hydrologic cycle of Sermeq 
Avannarleq. This tidewater glacier, which calves into a side 
arm of Jakobshavn Fiord, is located downstream of CU/ETH 
('Swiss') Camp. Its ice-dynamic flowline was determined by 



Fig. 2. Schematic of the model variables and coordinate system at a 
given node. Variable notation table is provided in the Appendix. 

following the path of steepest ice surface slope both up- and 
downstream of JAR2 automatic weather station (AWS; 
located at 69.42° N, 50.08° W in 2008) using a digital 
elevation model (DEM) of the Greenland ice sheet, with 
625 m horizontal grid spacing, derived from satellite altim- 
etry and enhanced by photoclinometry (Scambos and 
Haran, 2002). The resultant 530 km flowline runs up-glacier 
from the tidewater terminus of Sermeq Avannarleq (km 0 at 
69.37° N, 50.28° W), through JAR2 AWS (km 14), within 
2 km of Swiss Camp (km 46 at 69.56° N, 49.34° W in 2008), 
to the main ice divide of the Greenland ice sheet (km 530 at 
71 .54° N, 37.81 ° W; Fig. 1 ). The tidewater tongue of Sermeq 
Avannarleq was stable between 1955 and 1985 (Thomsen 
and others, 1988). Field observations suggest that since 
~2000 the tidewater tongue, currently ~2 km long, has 
retreated ~2 km and the ice surface of the terminal ~5 km of 
the glacier has become increasingly crevassed. 

2. METHODS 

2.1. Hydrology model 

We apply a 1-D (depth-integrated) hydrology model to the 
terminal 60 km of the Sermeq Avannarleq flowline to 
compute subglacial hydraulic head. We conceptualize the 
glacier hydrologic system as having perfect connectivity 
between the englacial and subglacial hydrologic systems 
(i.e. small hydraulic resistance). The assumption of hydraulic 
equilibrium between the englacial and subglacial systems is 
reasonable for capturing the seasonal timescale behavior of 
the glacier hydrologic system. In this conceptual model the 
subglacial hydraulic head is equivalent to the local englacial 
water-table elevation, h e . We therefore use the terms 
'subglacial hydraulic head' and 'englacial water-table 
elevation' interchangeably. Conduits, which operate at the 
ice/bed interface, whose geometry evolves through time, 
control the horizontal water discharge, Q, within the glacier 
hydrologic system (Fig. 2). Thus, h e varies in time and space 
due to variable conduit inflow and outflow as well as 
dynamic changes in conduit storage. As this model only has 
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one horizontal layer (or 'component') representing both the 
englacial and subglacial systems, it can be viewed as a 
simplification of more advanced multicomponent models 
that parameterize the supraglacial, englacial, subglacial and 
groundwater hydrology components independently (Flowers 
and Clarke, 2002; Kessler and Anderson, 2004). As a 
consequence of not explicitly representing the supraglacial 
hydrologic system and the hydraulic retention time therein 
(i.e. firn; Fountain and Walder, 1998), surface ablation is 
routed instantly to the top of the englacial water table. By 
enforcing water conservation, the rate of change in hydraulic 
head (or englacial water-table elevation) at a given node 
along the flowline may be calculated from four terms: 
(1) external meltwater input, / (via both surface and basal 
ablation per unit area), multiplied by the flowband cross 
width, w, (2) internal meltwater generation due to viscous 
melt within the conduit, m/p w , (3) horizontal divergence of 
conduit discharge, dQ/dx, and (4) changes in conduit 
storage volume (per unit length along the flowline) through 
time, dSJdt 


. dh e tin dQ dS c 

/iv + 

at p w ox at 


o) 


where 0 is the bulk ice porosity at a given node. With w= 1 , 
Equation (1) describes the flow per unit width along the 
flowline. Assuming the ice behaves like a fractured rock- 
type aquifer rather than a porous medium-type aquifer, this 
bulk ice porosity (or 'macro-porosity') is taken to represent 
the fractional volume occupied by fully connected surface- 
to-bed discrete water-storage elements (i.e. the lumped 
fractional volumes of surface and basal crevasses, conduits, 
moulins, etc.; Flowers and Clarke, 2002). The geometry of 
these discrete water-storage elements remains unspecified 
(Kessler and Anderson, 2004). Observed values of bulk ice 
porosity generally range between 0.004 and 0.013 in alpine 
glaciers (Fountain and Walder, 1 998). Similar to Kessler and 
Anderson (2004), we take bulk ice porosity as 0.01 and 
assume it is constant in space and time. We also run 
simulations with 0=0.005 and 0.015, the water content 
range constrained by ice temperature modeling of nearby 
Jakobshavn Isbrae (Luthi and others, 2002). 

In the conceptual framework of this model, the total 
water volume stored per unit length along the flowline at a 
given node within the glacier hydrologic system, 5, is the 
sum of two terms: 


S = S e + S c , (2) 

where S e is the englacial storage volume per unit length and 
5 C is the subglacial conduit storage volume per unit length. 
S e = cj)wH e , where F/ e is englacial water-column thickness 
(h e -h b ). Assuming a semicircular conduit geometry, conduit 
storage volume per unit length (or total cross-sectional area) 
can be expressed as a function of conduit radius, r: 

S c = -^-(n c w), (3) 

where n c is the number of conduits per unit width in the 
cross-flow (y) direction (defined in section 2.1.2). 


2.1 A . External meltwater input 

At each node along the flowline, external meltwater input, /, 
occurs as the sum of both surface ablation rate, a s , and basal 
ablation rate, a b , according to 

I — a s F + a b , (4) 


where F is the fraction of surface meltwater that reaches the 
englacial water table. This external meltwater input term 
does not accommodate episodic supraglacial lake drainage 
events, which can deliver tremendous amounts of water to a 
single location within the glacier hydrologic system (Das 
and others, 2008). Rather, our interest is in the seasonal 
timescale response to the annual surface ablation forcing. In 
the ablation zone, where no firn layer exists, almost all 
surface meltwater can be expected to reach the glacier 
hydrologic system. In the accumulation zone, however, 
meltwater must percolate vertically through the snowpack 
and firn. As a result of temporary storage or refreezing during 
this percolation, only a fraction of surface meltwater 
production reaches the glacier hydrologic system (Pfeffer 
and others, 1991; Fountain and Walder, 1998; Janssens and 
Huybrechts, 2000). We approximate the englacial entry 
fraction based on the ratio of annual surface accumulation, 
c s , to annual surface ablation, a s . This formulation follows 
Pfeffer and others (1991), who suggest that a given fraction 
of the annual accumulation must melt and saturate before 
runoff occurs, i.e. 

F=1-F r (VV) limit: F > 0. (5) 

V Iasi/ 

Here F r is the fraction of surface ablation retained in the firn 
where annual surface ablation and accumulation are equal. 
We take F r as 0.5, which implies that water enters the glacier 
hydrologic system upstream of the equilibrium line. The 
along-flowline profile of F varies relatively little when F r is in 
the range 0.3-0. 7 (Fig. 3). 

We interpolate along-flowline annual surface accumu- 
lation, c s , from a previously compiled dataset (Burgess and 
others, 2010). Annual surface ablation, a s , is taken as a 
function of elevation, based on previous in situ observations 
(Fausto and others, 2009), i.e. 

a s = ^h s -hf)-af, (6) 

where 7 is the present-day ablation gradient (Aa s /A/? s , taken 
as 0.00372; Fausto and others, 2009) and hf a is the 
equilibrium-line altitude (ELA, taken as 1125 m) and af a is 
the annual surface ablation at the ELA (taken as 0.4 m). 
Annual surface ablation is distributed through time to yield a 
surface ablation rate, a s , using a sine function to represent 
the melt-season solar insolation history (cf. Pimentel and 
Flowers, 2010): 

s > = 2 <7) 

where j is a given Julian date (JD), j 0 is the JD of melt onset 
(estimated as y 0 = 0.01 83/? s + 11 4), and D is the duration of 
the melt season (the JD of melt cessation is similarly 
estimated as ; c = -0.01 83F s + 248). While these idealized 
dates of melt onset and cessation lie within the observed 
range, this function likely underestimates the date of peak 
melt, which usually occurs closer to the end of the melt 
season (Fig. 4). 

The rate of basal ablation, a b , is calculated as 

• = (Qd + Qf + Qg) ( o\ 

b (Pi9 

where Qd and Qf are the heat input due to deformational 
(strain) and frictional heating, respectively, Q g is the 
geothermal flux, p\ is the density of ice (91 7 kg m -3 ) and L 
is the latent heat of fusion (333 550 J kg -1 ). We assume a 
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Fig. 3. (a) Observed annual surface accumulation, c s (Burgess and others, 2010), and ablation, a s (Fausto and others, 2009), versus 
distance upstream, (b) Englacial entry fraction, F, versus distance upstream with variable values of retention fraction, F r . Vertical dashed lines 
denote the locations of JAR2, Swiss Camp and the equilibrium line, (c) Englacial entry fraction (with F r = 0.5) over the observed range of c s 
and a s values. 


geothermal flux of 57mWm“ 2 along the entire flowline 
(Sclater and others, 1980). We also assume that deforma- 
tional heating is concentrated at the ice/bed interface where 
it contributes to basal ablation (EHooke, 2005). The lack of 
vertical resolution in the single-component hydrology model 
also necessitates that internally produced meltwater is routed 
instantaneously to the bed (Flowers and others, 2005). 
Deformational and frictional heating are calculated as 


Qd + Qf — r b( u + u b )/ (9) 

where r b is taken as the basal gravitational driving stress, u 
is depth-averaged deformational ice velocity and u b is basal 
sliding velocity (EHooke, 2005). Both r b and u are calcu- 
lated according to the shallow-ice approximation using the 
flow-law parameter for ice defined in section 2.1.4. As this 
glacier hydrology model is not currently coupled to an ice 
flow model that incorporates transient basal sliding vel- 
ocities, we simply parameterize u b as 1 5 m a -1 throughout 
the ablation zone. Based on observations, this basal sliding 
velocity can be regarded as characteristic of the Sermeq 
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Fig. 4. The modeled time-space distribution of surface ablation 
rate, |a s |. Vertical dashed lines denote the locations of JAR2, Swiss 
Camp and the equilibrium line. 


Avannarleq ablation zone (see Colgan and others, in 
preparation). 

First-order calculations suggest that Qf and Q g are small 
compared with Qd, which makes Qd the main control on 
basal ablation. 

With a mean value of 3.3 cm (std dev. 1.6cm) in the 
ablation zone, annual basal ablation is a very small fraction 
of annual surface ablation. Beneath the km 7 icefall, 
however, steep surface slopes produce a strong deforma- 
tional heat flux, which results in annual basal ablation of 
9.3 cm. In comparison, annual surface ablation at the icefall 
is ~3.5m. These annual basal ablation values are small 
compared with those estimated for Jakobshavn Isbrae, where 
deformational heating alone is capable of producing >0.5 m 
of meltwater annually (Truffer and Echelmeyer, 2003). 
Unlike surface ablation, basal ablation is distributed evenly 
throughout the year. We neglect submarine basal ablation at 
the ice/ocean interface, as hydraulic head is prescribed as 
sea level within the ice tongue. 


2. 1.2. Conduit discharge 

Horizontal divergence of conduit discharge, dQ/dx, is the 
along-flowline gradient in conduit discharge, Q. For sub- 
glacial conduits, where Reynolds numbers are expected to 
be large, we calculate discharge based on the Darcy- 
Weisbach equations for turbulent flow in conduits (e.g. 
Rothlisberger, 1972; Nye, 1976; Spring and Hutter, 1982; 
Flowers and others, 2004; Pimentel and Flowers, 2010). As 
we assume perfect connectivity between the subglacial and 
englacial systems and that conduits are located at the ice/ 
bed interface, conduit discharge is dependent on the local 
hydraulic-head gradient, dhjdx\ 




dh e 
~d x 


(n c w) sign 



( 10 ) 


where g is gravitational acceleration (9.81 ms -2 ), f is a 
friction factor, D h is the effective hydraulic diameter, r is the 
conduit radius, n c is the number of conduits per meter in the 
cross-flow direction and w is the flowband cross width. We 
assume a friction factor of 0.05, which is appropriate for 
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Fig. 5. (a) Total annual external meltwater input per unit width 
entering the glacier hydrologic system upstream of a given position 
(f I • dx). (b) Maximum conduit radius, r max . (c) Conduit spacing in 
the across-flow direction, n c . (d) Total conduit storage volume per 
unit length, S c . Line color varies with glacier hydrology length 
scale, a. Vertical dashed lines denote the locations of JAR2, Swiss 
Camp and the equilibrium line. 

turbulent flow in rough-walled pipes, where the amplitude 
of the wall roughness elements is ~0.02 times the pipe 
diameter (Moody, 1 944). To assess model sensitivity, we also 
present a simulation in which f is taken as 0.01. 

We approximate the observed arborescent nature of 
subglacial conduit systems (Fountain and Walder, 1998) 
within the framework of our flowline model by making the 
number of conduits per unit width, n c , in the cross-flow (y) 
direction an exponential function of distance upstream from 
the terminus, x: 

n c = n' er exp0), (11) 

where r£ er is the number of conduits per unit width at the 
terminus (taken as 0.005 or 200 m between conduits) and a is 
a glacier hydrology length scale controlling the rate of 
increase in number of conduits per unit width (or the rate of 
decrease in spacing between conduits) with distance up- 
stream. At any flowline position, the subglacial hydrologic 
system is represented as a set of non-interacting parallel 
conduits, with a conduit spacing given by Equation (1 1 ). 

As disequilibrium between the forces promoting conduit 
opening and closure is permitted in the model (discussed in 
section 2.1.4), we must constrain conduit radii between 
reasonable maximum and minimum dimensions to prevent 
implausible feedbacks (i.e. permanent closure and unre- 
stricted growth). The need to impose these boundaries may 
be attributed to the intrinsic limitations of a simplified 1-D 
hydrologic model. The evolving geometry of the subglacial 
hydrologic system can only be represented accurately with a 
higher-dimensional model that allows for flow divergence 
(convergence) away from (towards) highly conductive con- 
duits. A model that incorporates these processes can in 
principle simulate conduit geometry without prescribed 



Fig. 6. Schematic of the arborescent approximation of the conduit 
network: number of conduits (rounded to nearest integer) and 
maximum conduit diameter, 2 r (5 times exaggerated), versus 
distance upstream when a = 20 km. Horizontal dashed lines denote 
the locations of JAR2, Swiss Camp and the equilibrium line. 


upper and lower bounds for conduit radii. Higher-dimen- 
sional modelling would also benefit from a more precise 
characterization of the bedrock geometry throughout our 
study area. 

Similar to conduit spacing, we also parameterize max- 
imum conduit radius, r max , as an exponential function of 
distance upstream from the terminus: 

''max = C.x eX p(“^)' O 2 ) 

where r^ x is the maximum permissible conduit radius at the 
terminus (arbitrarily taken as 2 m) and a is the same glacier 
hydrology length scale as in Equation (11), which controls 
the rate of decrease in maximum conduit radius with 
distance upstream (r min is taken as 10% of r max ). 

On the assumption that the total conduit storage volume 
per unit length, 5 C , at a given position along the flowline 
reflects the integrated annual external meltwater input per 
unit width that enters the glacier hydrologic system up- 
stream of that position (f I dx), we can approximate the 
glacier hydrology length scale, a, as the length scale of the 
along-flowline f I • dx profile. A qualitative assessment 
suggests that a = 20 km renders an along-flowline S c profile 
that approximates the f / • dx profile better than o = 10 or 
30 km (Fig. 5). Thus, we take a = 20 km as the glacier 
hydrology length scale. Together, Equations (11) and (12) 
capture the essence of an arborescent conduit network, 
namely the decrease in both conduit spacing and conduit 
radius with distance upstream (Fig. 6). There are, however, 
an infinite number of combinations of conduit spacing and 
radius values that may yield a given discharge, and neither 
parameter is tightly constrained by field observations. Thus, 
our maximum and minimum conduit radii may not 
necessarily reflect actual limits found in nature; they merely 
reflect limits that, together with the conduit spacing relation 
we have imposed, lead to reasonable values for the overall 
transmissivity of the conduit system. 
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Fig. 7. (a) Modeled (black; see Colgan and others, in preparation) 
and observed (grey; Scambos and Haran, 2002) ice surface 
elevation, h s , and observed bedrock topography, h b (Bamber and 
others, 2001; Plummer and others, 2008). (b) Estimated basal ice 
temperature, 7j. Vertical dashed lines denote the locations of JAR2, 
Swiss Camp and the equilibrium line. 


2.1.3. Meltwater generation in conduits 
At each node along the flowline, meltwater is generated 
within the conduits as a consequence of the dissipation of 
the heat generated by viscous friction in the water; this heat 
is instantaneously conducted to conduit walls (Nye, 1976). 
The rate at which water is produced by internal melt, rn, in 
units of mass per length per time at a given node, is a 
function of conduit discharge, Q, local hydraulic-head 
gradient, dhjdx, and the latent heat of fusion of ice, L: 




dh e 

dx 


(13) 


2.1.4. Conduit volume 

As the conduits grow and shrink at each node along the 
flowline, the local conduit storage volume per unit length, 5 C 
(or total cross-sectional area), changes. Conduit storage 
volume grows through two processes: (1) the volume per 
unit length created by the ice removed through internal 
meltwater generation (m/p\, where m is calculated according 
to Equation (13)), and (2) deformational opening, which 
occurs when water pressure, P w (defined as p w gH e ), exceeds 
ice pressure, P\ (defined as p\gH\). Conversely, conduit 
storage volume only decreases through deformational clos- 
ure when P\ > /V Thus, our model assumes that deforma- 
tional opening of conduits can occur in response to high 
subglacial water pressures (Nye, 1976; Ng, 2000; Clarke, 
2003). As we are assuming perfect connectivity between the 
subglacial and englacial systems, the basal water pressure 
both inside and outside the conduits may be regarded as 
equivalent. This implicitly discounts the possibility of 
leakage' of high water pressures from the conduit system 
to adjacent regions (Bartholomaus and others, 2008). As 
basal ice temperature in the ablation zone is likely at the 
pressure-melting point (Phillips and others, 201 0), we do not 
consider conduit closure through refreezing of water within 
the conduits. The processes that govern the transient rate of 
change of conduit storage volume per unit length, OSJdt, 
can be expressed as (Nye, 1976) 


dSc 

dt 



(■ n c w) ( ^— n Pw ^ sign(Pi - P w ), 


(14) 


where A is the temperature-dependent flow-law parameter 


(EHuybrechts and others, 1991) and n is the Glen's law 
exponent (n = 3). We enhance the temperature-dependent 
flow- 1 aw parameter by a factor of three to account for softer 
Wisconsinan basal ice (Reeh, 1985; Paterson, 1991), which 
comprises the bottom ~20% of the ice column along the 
flowline (Huybrechts, 1994). 

As recent thermodynamic modeling suggests that basal 
ice throughout the ablation zone of the Sermeq Avannarleq 
flowline is at the pressure-melting point (Phillips and others, 
2010), we calculate basal ice temperature, 7j, downstream 
of the equilibrium line as 

T\ = 273.1 5 — /3H\, (15) 

where [3 is the change in melting point with ice thickness 
(taken as 8.7 x 10 _4 Km _1 ; Paterson, 1994). In the accumu- 
lation zone, where downward vertical ice velocities are 
expected to produce cold basal ice temperatures, we simply 
prescribe basal ice temperature as 268 K. This is the basal ice 
temperature predicted just upstream of the equilibrium line 
by recent thermodynamic modeling (Phillips and others, 
unpublished information). We linearly transition from pres- 
sure-melting point basal ice temperatures in the ablation 
zone to the cold basal ice temperature of the accumulation 
zone between the equilibrium line and 5 km upstream of the 
equilibrium line. The resultant basal ice temperature profile 
is consistent with the notion that meltwater percolation 
rapidly warms the ice below the equilibrium line and 
captures the essence of Sermeq Avannarleq basal ice 
temperature profiles produced by Phillips and others 
(unpublished information; Fig. 7). 

2.2. Datasets and boundary conditions 

This hydrology model employs several observational data- 
sets. Annual surface accumulation and annual surface 
ablation profiles were obtained from previously compiled 
datasets (Fausto and others, 2009; Burgess and others, 201 0). 
Both these ice-sheet-scale surface climatology datasets use 
models to interpolate between in situ observations, in- 
cluding the four Greenland Climate Network (GC-Net) 
AWSs along the Sermeq Avannarleq flowline (Steffen and 
Box, 2001). To define ice geometry, we use a modeled ice 
surface elevation (h s ) profile (see Colgan and others, in 
preparation), which compares well with the profile inter- 
polated from a DEM of the Greenland ice sheet derived from 
satellite altimetry enhanced by photoclinometry (Scambos 
and Haran, 2002; Fig. 7). A coarse-resolution bedrock 
elevation (h b ) dataset is available along the entire flowline 
(Bamber and others, 2001). While this dataset has nominal 
5 km horizontal resolution, actual resolution is dependent 
on the spacing of airborne flight-lines, and thus varies 
spatially and can exceed 5 km. A finer-resolution dataset 
(nominally 750 m) is available for the terminal 40 km of the 
flowline (Plummer and others, 2008). 

Equipotential englacial hydraulic-head gradients are 
theoretically ~1 1 times more sensitive to ice surface slopes 
than bedrock slopes (Shreve, 1972). Thus, first-order ice 
surface slope (i.e. regional or ~10 ice thicknesses) may be 
expected to govern the first-order geometry of subglacial 
flow direction. For example, the water beneath Swiss Camp 
can be expected to travel ~30 km to the ice-sheet margin on 
an azimuth of ~279° (i.e. following regional ice surface 
slope; dashed line in Fig. 1). By forcing subglacial water to 
flow along the ice-dynamic flowline to the margin (i.e. 
following steepest local or ~1 ice thickness surface slope), 
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Fig. 8. Total flowline prescribed external water input, f Iw • dx 
(m 3 a -1 ), and total flowline modeled water storage, fS- dx, for 
various values of bulk ice porosity, 0, during 20 year simulations. 

subglacial water is routed ~48 km on an azimuth of ~225°. 
This artificially reduces the hydraulic-head gradient, dhjdx, 
by artificially increasing dx. To compensate for this we use 
simple trigonometry to project the dynamic flowline on the 
hydrologic coordinates, correcting dx values on a node-by- 
node basis, and thereby maintaining realistic dhjdx values 
along the length of the flowline. 

We apply first-type (specified head) Dirichlet boundary 
conditions to hydraulic head, h e . From kmO to 2, the 
tidewater terminus, h e is prescribed as sea level (i.e. 
h e = 0m). We use the cold-to-warm basal ice temperature 
transition at the equilibrium line (km 51 ) to approximate the 
upstream limit of the glacier hydrologic system. We impose a 
no-flux boundary that extends upwards from the bedrock 
beneath the equilibrium line to the ice surface at 1 1 times the 
local surface slope. This is the steepest theoretical gradient 
that water traveling through the ice sheet could be expected 
to take in order to reach the bed at the cold-to-warm basal ice 
temperature transition beneath the equilibrium line (Shreve, 
1972). Thus, upstream of the equilibrium line, there is a 
portion of the glacier hydrologic system that behaves as a 
'perched' aquifer (i.e. underlain by cold ice). In this ~4km 
portion of the flowline, we assume that ice temperature is 
sufficiently heterogeneous to allow the persistence of 
englacial conduits (Catania and Neumann, 2010). We treat 
the discharge of these conduits in an identical manner to the 
conduits downstream of the equilibrium line, except we 
suspend the processes that govern the transient rate of change 
of conduit storage volume, dSJdt , and prescribe conduit 
radii as constant (5 cm). We assume that any water entering 
the ice upstream of the intersection of the perched aquifer 
and the ice surface refreezes, although we do not perform a 
full enthalpy solution to incorporate this process as a 
transient flux. We also assume that basal ablation is zero 
beneath the cold basal ice in the accumulation zone. The 
glacier hydrologic system is therefore fully transient over the 
49 km between the equilibrium line (km 51) and the tide- 
water terminus (km 2). 

The differential equations describing transient hydraulic 
head, dhjdt, were discretized in space using first-order finite- 
volume methods (dx= 500 m) with hydraulic head, h e , at cell 
centers and fluxes, Q, at cell edges. The semi-discrete set of 
ordinary differential equations coupled at the computational 
nodes was then solved using 'ode15s', the stiff differential 
equation solver in MATLAB R2008b, with a 1 day time-step. 



x (km) 


Animation 1. Modeled hydrology along the terminal 60 km of the 
Sermeq Avannarleq flowline, (a) External meltwater input, I; both 
surface (red; a s F) and basal ablation (magenta; ab; 30 x exagger- 
ated). (b) Ice surface elevation (h s ; white), bedrock elevation 
(/?b; brown) and hydraulic head (or englacial water-table elevation; 
h e ; blue). Dashed green line represents the hydraulic head equal to 
flotation, (c) Conduit radius, r. (d) Rate of change in hydraulic head, 
dh e /dh t . Vertical dashed lines identify the positions of JAR2, Swiss 
Camp and the equilibrium line. Full movie available at www. 
igsoc.org/hyperl ink/1 0J1 54_Animation1 .mov. 


3. RESULTS 

The 1-D hydrology model achieves a quasi -steady-state 
annual cycle in hydraulic head after years of spin-up 
(Fig. 8; Animation 1). Generally, modeled flotation fraction, 
P w /Pj, increases with distance upstream to a maximum at the 
equilibrium line (Fig. 9). Upstream of the equilibrium line, 
where basal ice temperatures are assumed to transition from 
warm to cold, the englacial hydrologic system likely 
behaves as a perched aquifer where it becomes underlain 
by cold ice. In this region, the high englacial water-table 
elevation should not be interpreted as a high flotation 
fraction, because this water pressure is not exerted at the 
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Di sta nee u pst rea m ( km ) 

Fig. 9. Modeled time-space distribution of flotation fraction, PJP\. 
The white contour denotes a flotation fraction of 1 . Vertical dashed 
lines denote the locations of JAR2, Swiss Camp and the 
equilibrium line. 
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Fig. 10 . Modeled time-space distribution of the rate of change of hydraulic head (or englacial water-table elevation; dhjdt ) when conduit 
friction factor, f, equals 0.05 (a) and 0.01 (b). Color bar saturates below -3.6 m d _1 . Vertical dashed lines denote the locations of JAR2, Swiss 
Camp and the equilibrium line. 



bed. In fact, in this region, we assume that no subglacial 
water is present and hence subglacial water pressure should 
be zero. The km tidewater tongue remains at flotation 
year-round (i.e. P w /P\ = ^). In response to meltwater input 
from surface ablation, flotation fraction exhibits a broad 
summer peak throughout the ablation zone. This peak 
occurs progressively later with distance upstream, due to the 
upstream migration of surface ablation inputs. The modeled 
hydraulic head fluctuates relatively close to flotation 
throughout the year; flotation is achieved for a brief period 
between km 1 6 and 33. 

The time-space pattern of the rate of change in hydraulic 
head, dhjdt (a proxy for the rate of change in glacier water 
storage, d5/df), may be conceptualized as the derivative of 
the time-space pattern of flotation fraction (when the 
condition dS c /dt < dS e /dt is satisfied; Fig. 10). The entire 
flowline experiences positive rates of glacier water storage 
(or increasing hydraulic head) at the beginning of the melt 
season (as high as 3.6 md -1 ), and negative rates of glacier 
water storage (or decreasing hydraulic head) at the end of 
the melt season (as low as -10.2 md -1 ). Generally, the total 
water volume stored along the flowline decreases (or drains) 
more slowly than it increases (or fills; Fig. 8). The magnitude 
of the modeled annual dhjdt cycle decreases with distance 
upstream, reaching approximately zero in the vicinity of the 
equilibrium line (taken as the upstream boundary of 
temperate basal ice and the englacial hydrologic system). 
The transition from positive to negative rates of storage also 
progressively lags with distance upstream. Unlike the 
positive rates of glacier water storage, which exhibit a 
smooth time-space distribution, the negative rates of glacier 
water storage have a more patchy time-space distribution, 
with certain dates and flowline locations experiencing 
anomalously fast or slow negative rates of glacier water 
storage. This is a consequence of the evolution of the 
subglacial conduit system. 

Modeled subglacial conduit storage exhibits sharp along- 
flowline boundaries between regions of closed and opened 
conduits (i.e. km 7, 16, 25, 42 and 50; Animation 1). This 
suggests two types of subglacial environments exist: (1 ) 'year- 
round' (i.e. km 16-25 and 42-50) and (2) 'seasonally' (i.e. 
km 7-1 6 and 25-42) open subglacial conduit storage. 
Following the onset of surface ablation, conduits begin to 
grow upstream from the tidewater terminus (where conduits 


are open year-round due to a lack of deformational closure 
pressure in the floating tidewater tongue). When a stretch of 
seasonally open subglacial conduit storage connects to an 
upstream stretch of year-round open subglacial conduit 
storage, there is a temporary anomalously large negative rate 
of water storage (i.e. decrease in local hydraulic head) as the 
portion of the flowline underlain by the year-round conduits 
is drained relatively quickly. Subsequently, the upstream 
stretch of seasonally open subglacial conduit storage begins 
to open (due to an increase in local hydraulic-head gradient) 
and the negative rate of water storage returns to a smaller 
magnitude. The opening of seasonal conduits that serve to 
connect stretches of year-round conduit storage are respon- 
sible for the patchy time-space distribution of negative rates 
of glacier water storage. 

4. DISCUSSION 
4.1. Hydrology model 

Based on the absolute value of flotation fraction, as well as 
qualitative features in the annual hydraulic-head cycle (i.e. a 
positive phase followed by a negative phase with no 
significant change during the winter), we conclude that 
the relatively simple one-component 1-D hydrology model 
produces a reasonable annual hydrologic cycle along the 
terminal 60 km of the Sermeq Avannarleq flowline. The 
modeled hydraulic head oscillates close to flotation 
throughout the ablation zone of the Sermeq Avannarleq 
flowline. The suggestion of a relatively high englacial water 
table is consistent with previous in situ observations of 
flotation fractions between 0.79 and 1.05 in the Paakitsoq 
area (69.5° N, 50.0° W; Thomsen and Olesen, 1991) and 
between 0.95 and 1.00 in Jakobshavn Isbrae (69.2° N, 
48.7° W; Iken and others, 1993; Luth i and others, 2002). The 
high modeled flotation fraction along the flowline is due to 
the dependence of conduit conductivity and discharge on 
the relative difference between water and ice pressures. 
Sufficiently high hydraulic head is needed to overcome 
basal ice pressure and open the conduits. Once open, 
conduits draw down hydraulic head until reaching a 
threshold below which hydraulic head is no longer 
sufficient to counteract basal ice pressure and keep the 
conduits open. The conduits subsequently close, achieving 
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Fig. 11. Mean annual flotation fraction, P w /Pj, versus distance 
upstream for various values of bulk ice porosity, 0. Dashed curves 
represent annual maximum and minimum values. Vertical 
dashed lines denote the locations of JAR2, Swiss Camp and the 
equilibrium line. 


their minimum prescribed radius, r min , leaving large 
amounts of residual water stored in the glacier. 

Model runs in which bulk ice porosity, 0, is taken as 
0.005 and 0.015 return very similar along-flowline profiles 
of mean annual flotation fraction (Fig. 1 1). Thus, a relatively 
high hydraulic head is likely an inherent feature of a glacier 
hydrologic system drained solely by conduits. The along- 
flowline profile of maximum flotation fraction is, however, 
higher (lower) in the lower (higher) bulk ice porosity 
simulation. Thus, as currently parameterized, the magnitude 
of dhjdt is a function of bulk ice porosity. When 0 = 0.005, 
flotation fraction values are unrealistically high in the 
vicinity of km 25 (~1 .2 which is conducive to artesian 
springs). Bulk ice porosity values as low as 0.0005 have 
been postulated for alpine glaciers (Humphrey and others, 
1986). Our sensitivity analysis suggests that increased water 
storage at the bed must be invoked to produce realistic 
flotation fraction values when 0 <c 0.01 . At present, our 
single-component model does not explicitly incorporate the 
water storage volume created at the ice/bed interface by 
vertical uplift of the ice sheet (cf. Pimentel and Flowers, 
2010). Due to this limitation a portion of the bulk ice 
porosity can be assumed to represent the storage volume of 
discrete voids at the ice/bed interface. In reality, the storage 
volume at the ice/bed interface is likely a function of 
hydraulic head, rather than constant in space and time as we 
have implicitly assumed by using constant porosity. 

The model inference of year-round subglacial conduit 
storage along certain portions of the flowline (km 1 6-25 and 
42-50) is supported by limited observations of the persist- 
ence of some englacial hydrology features through the 
winter in the Sermeq Avannarleq ablation zone (Catania and 
Neumann, 201 0). The sharp model transitions between year- 
round and seasonally active conduit storage appear to be 
due to slight changes in the along-flowline gradient in local 
hydraulic head, dhjdx ; subglacial conduit storage volume 
per unit length, 5 C , is proportional to the magnitude of {dhj 
ax) 1 ' 5 . When a bedrock topography that increases mono- 
tonically upstream is imposed, the pockets of year-round 
subglacial conduit storage form in the same location as 
when the observed bedrock topography is imposed (Anima- 
tion 2). As short stretches of year-round subglacial conduit 
storage appear to form downstream of stretches of relatively 
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Animation 2. Same as Animation 1 except with a modified bedrock 
topography that increases monotonically upstream. Full movie 
available at www.igsoc.org/hyperlink/10J154_Animation2.mov. 


steep ice surface topography, we speculate that the along- 
flowline gradient in ice surface elevation, dhjdx, influences 
the along-flowline gradient in local hydraulic head, dhjdx, 
and hence the location of year-round subglacial conduit 
storage. When an ice surface topography with sinusoidal 
undulations (25 m amplitude and 10 km wavelength) is 
imposed, the pockets of year-round subglacial conduit 
storage tend to form beneath portions of the flowline with 
the steepest ice surface slopes (Animation 3). Along-flowline 
gradients in flotation fraction, P w /P\, as a result of ice surface 
topography may also have an influence in the location of 
year-round subglacial conduit storage. Maintaining year- 
round subglacial conduit storage through deformational 
opening, however, could only be expected along portions of 
the flowline where water pressure exceeds ice pressure year- 
round. The time-space distribution of modeled flotation 
fraction suggests that this only occurs in the tidewater 
tongue (km 0-2; Fig. 9). The model inference that portions of 
the subglacial conduit storage system may overwinter at the 
base of the Greenland ice sheet to be reactivated the 
following melt season differs slightly from previous alpine 
glacier work, which suggests a new channelized subglacial 
hydrologic system migrates up-glacier from the terminus 
each melt season (Hubbard and Nienow, 1997). 

During the period of peak decreasing hydraulic head at 
Swiss Camp (i.e. JD 300 or 325, depending on /"parameter- 
ization; Fig. 10) the conduits draw down hydraulic head 
with the efficiency of a semicircular conduit 40cm in 
diameter spaced every 20 m and operating at basal ice and 
water pressure (i.e. n c = 0.050 m and r= 0.20 m). Similarly, 
during the period of peak decreasing hydraulic head at 
JAR2 HD 180) the conduit system is parameterized to 
discharge water with the efficiency of a semicircular conduit 
1.98 m in diameter spaced every 100 m operating at basal 
ice and water pressure (n c = 0.010 m and r= 0.99 m). 
Although our parameterized conduit system appears to be 
reasonable, a lack of field observations makes this config- 
uration speculative. The choice of conduit friction factor, f, 
which is a measure of resistance (discharge is inversely 
related to f), appears to influence the timing, but not 
magnitude, of negative dhjdt values (Fig. 10). A smaller 
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Animation 3. Same as Animation 2 except with a sinusoidal pattern 
(25 m amplitude and 10km wavelength) imposed on ice surface 
topography and basal ablation actively disabled. Full movie 
available at www.igsoc.org/hyperlink/10J154_Animation3.mov. 


value of f propagates the perturbation in hydraulic-head 
gradient responsible for conduit opening upstream more 
quickly than a larger value of f. 

Assuming a bulk ice porosity of 0.01, and that the 
modeled glacier hydrologic cycle is in quasi-steady state, 
the mean residence time, f res , of water in the flowline can 
be calculated as ~2. 2 years, when w=1 m, according to 


JS - d * , 

J Iw - dx' 


(16) 


where f S • dx is the total flowline modeled water storage 
(annual mean ^45 000 m 3 ) and f Iw • dx is the total 
flowline prescribed external water input (annual mean 
^20 000 m 3 a -1 ; Fig. 8). The integral notation (f -dx) refers 
to integration of continuous values over the entire length of 
flowline (i.e. all x values). Assuming bulk ice porosities of 
0.005 and 0.015, the mean residence times are ~1.1 and 
3.3 years, respectively. This relatively short residence time 
suggests that the glacier hydrologic system can be expected 
to respond to external meltwater forcings on a relatively 
short timescale (i.e. years rather than decades). 


4.2. Relevance to basal sliding 

In section 1 we reviewed the notion that variations in basal 
sliding velocity can be attributed to variations in the rate of 
change of glacier water storage (i.e. dS/dt) in alpine glaciers 
(Kamb and others, 1994; Fountain and Walder, 1998; 
Anderson and others, 2004; Bartholomaus and others, 
2008). This conceptual model of basal sliding has recently 
been extended to the Greenland ice sheet (Bartholomew 
and others, 2010). Below we present both remotely sensed 
and in situ velocity data from the Sermeq Avannarleq 
flowline that suggest that observed periods of enhanced 
basal sliding generally correspond to modeled periods of 
positive rates of change in glacier water storage (or 
increasing hydraulic head), while periods of reduced basal 
sliding conversely correspond to modeled periods of nega- 
tive rates of change in glacier water storage (or decreasing 
hydraulic head). 
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Fig. 1 2. Annual ice surface velocity cycle along the terminal 60 km of 
the Sermeq Avannarleq flowline (color bar saturates at 175 m a -1 ). 
Dotted black lines indicate individual InSAR velocity profiles. 
Vertical dashed lines denote the locations of JAR2, Swiss Camp 
and the equilibrium line. 


Given the proximity of Sermeq Avannarleq to Jakobshavn 
Isbrae, multiple 2005 and 2006 interferometric synthetic 
aperture radar (InSAR)-derived ice surface velocity profiles 
are available for the terminal 60km of the Sermeq 
Avannarleq flowline (Joughin and others, 2008). Linear 
interpolation between these profiles yields a time-space 
velocity plot (Fig. 12). While summer velocities are not 
available upstream of km 22, the region around JAR2 
exhibits a distinct annual velocity cycle. This cycle consists 
of a summer speed-up event, in which velocities exceed 
winter velocity, followed by a fall slowdown event, in which 
velocities decrease to lower than winter velocity. These 
observed velocity anomalies may be interpreted as periods 
of enhanced and suppressed basal sliding that qualitatively 
correspond to periods of modeled increasing and decreasing 
hydraulic head. 

FHigh-frequency differential GPS measurements have 
been acquired at Swiss Camp since June 1996, providing 
daily resolution of the annual ice surface velocity cycle for 
the last 12 years (Larson and others, 2001; Zwally and 
others, 2002; Fig. 13). In accordance with the InSAR 
observations, these in situ GPS observations reveal that the 
ice at Swiss Camp also experiences a summer speed-up 
event, in which ice velocities increase above winter 
velocities, followed by a fall slowdown event, in which ice 
velocities decrease below winter velocities. The summer 
speed-up event varies in peak magnitude and timing over 
the 12 year record, presumably due to variable surface 
climatology and meltwater forcing, as well as local 
supraglacial lake drainage events. While the model is not 
forced with an observed melt history, the modeled period of 
increasing hydraulic head (i.e. JD 150-200) generally 
matches the timing of the observed speed-up event. The 
observed slowdown event, however, follows the observed 
speed-up event in relatively quick succession. The modeled 
peak in decreasing hydraulic head lags the modeled peak in 
increasing hydraulic head by 125-150 days (depending on 
the value of f; Fig. 1 0). 

We speculate that this apparent delay in the negative dhj 
dt phase, which is due to a delay in opening the upstream 
conduits, is a consequence of 1-D hydrological modeling. A 
two-dimensional (2-D; xy) hydrological model would better 
capture changes in englacial water-table gradient at Swiss 
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Fig. 13. Observed GPS annual velocity cycle at Swiss Camp over the 1996-2008 period. 


Camp by allowing conduit discharge both parallel and 
perpendicular to the ice-dynamic flowline. Presently, by 
forcing water to flow along the entire ice-dynamic flowline 
to the terminus, the 1-D model essentially requires the 
perturbation in local hydraulic head gradient that initiates 
conduit opening to slowly migrate from the terminus 
upstream to Swiss Camp. In reality, the glacier hydrologic 
system near Swiss Camp is governed by both along- and 
across-flow hydraulic head gradients ( dhjdx and dhjdy 
respectively). A 2-D (xy) model would better represent the 
influence of bedrock topography on subglacial discharge, by 
allowing subglacial discharge to concentrate in topographic 
low points (Shreve, 1 972; Flowers and others, 2005). Due to 
this concentration of discharge, 2-D models are inherently 
more capable of developing larger and more efficient 
conduits than 1-D models. 

Based on the qualitative agreement between observations 
of a strong annual ice velocity cycle and the modeled 
annual glacier hydrologic cycle, we propose that the 
terminal 50 km of the Greenland ice sheet's Sermeq 
Avannarleq flowline behaves akin to an alpine glacier, 
whereby the summer speed-up event is caused by inefficient 
subglacial drainage and positive rates of change of glacier 
water storage (i.e. dS/df>0, or increasing hydraulic head), 
while the fall slowdown event reflects the establishment of 
efficient subglacial drainage and negative rates of change of 
glacier water storage (i.e. d5/df<0, or decreasing hydraulic 
head). A notable departure from the alpine analogy, 
however, is the model inference that subglacial conduit 
storage may be capable of overwintering beneath the 
Greenland ice sheet to be reactivated the following melt 
season, rather than migrating up-glacier from the terminus 
each melt season as observed in alpine glaciers (Hubbard 
and Nienow, 1997). The recovery of both basal sliding 
velocities and rates of change of hydraulic head, dhjdt , to 
winter values shortly after the fall slowdown, suggests that 
the glacier hydrologic system is capable of repressurizing to 
maintain winter sliding. Presumably, this occurs by: (1) a 
reduction in transmissivity due to conduit closure and (2) 
continued meltwater input, largely via basal ablation due to 
geothermal, frictional and deformation (strain) heat fluxes 
rather than internal meltwater generation, m/p w . 


5. SUMMARY REMARKS 

We developed a relatively simple one-component 1-D 
hydrology model to track glacier water storage and dis- 
charge through time. In this model, glacier water input is 
prescribed based on observed ablation rates, while glacier 
water output occurs through conduit discharge. Conduit 
discharge varies in response to the dynamic evolution of 
conduit radius. The hydrology model suggests that sub- 
glacial hydraulic head (or, equivalently in this model, 
englacial water-table elevation) annually oscillates relatively 
close to flotation, even reaching flotation for brief periods 
along certain stretches of the flowline. This is consistent with 
the few available borehole observations of englacial water- 
table elevation. Alternatively imposing idealized bedrock 
and ice surface topographies suggests that along-flowline 
gradients in ice surface elevation (or ice surface slope), 
rather than along-flowline gradients in bedrock elevation (or 
bedrock slope), control the locations of year-round sub- 
glacial conduit storage. The sharp transitions between 
subglacial environments that support year-round and sea- 
sonally open subglacial conduit storage are likely due to 
slight along-flowline changes in the local englacial 
hydraulic head gradient. A calculated mean glacier water 
residence time of ~2. 2 years implies that large amounts of 
water are stored in the glacier throughout the year. The 
glacier hydrologic system can therefore be expected to 
respond to external meltwater forcings (i.e. reorganize) on a 
relatively short timescale. A qualitative comparison between 
the observed annual ice velocity cycle and the modeled 
annual cycle of glacier water storage suggests that enhanced 
(suppressed) basal sliding occurs during periods of positive 
(negative) rates of glacier water storage. Thus, we speculate 
that the terminal 50 km of the Sermeq Avannarleq flowline 
likely experiences a basal sliding regime similar to that of an 
alpine glacier. Considering the inherent limitations of 1-D 
modeling, the timing of the modeled increasing and 
decreasing hydraulic-head events is also reasonable. Work 
is underway to produce a 2-D version of this hydrology 
model, which removes the need for a prescribed conduit 
geometry and allows a more refined characterization of the 
annual hydrologic cycle in the Sermeq Avannarleq region 
and other portions of the Greenland ice sheet. 
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APPENDIX 

Variable notation 

a Glacier hydrology length scale (m) 

(3 Change in melting point with ice thickness (Km -1 ) 
7 Ablation gradient 

0 Bulk ice porosity 

Pi Density of ice (kgrrT 3 ) 

p w Density of water (kg m -3 ) 


r\j Gravitational driving stress (Pa) 

A Flow law parameter (Pa -3 a -1 ) 

D Melt-season duration (days) 

Dh Effective hydraulic diameter (m) 

F Entry fraction 

F r Retention fraction 

H e Englacial water thickness (m) 

H\ Ice thickness (m) 

/ Rate of external meltwater input (ma _1 ) 

L Latent heat of fusion of ice (J kg -1 ) 

Pj Ice pressure (Pa) 

P w Water pressure (Pa) 

Q Conduit discharge (m 3 a -1 ) 

Qd Heating due to deformation (strain) (WrrT 2 ) 

Qf Heating due to friction (W rrT 2 ) 

Q g Geothermal flux (Wm -2 ) 

5 Total storage per unit length (m 2 ) 

5 C Subglacial conduit storage volume per unit length 

(m 2 ) 

5 e Englacial storage volume per unit length (m 2 ) 

T\ Basal ice temperature (K) 

ab Basal ablation rate (ma -1 ) 

a s Annual surface ablation (m) 

a s Surface ablation rate (ma -1 ) 

c s Annual surface accumulation (m) 

c s Surface accumulation rate (ma -1 ) 

f Conduit friction factor 

g Gravitational acceleration (ms -2 ) 

h b Bedrock elevation (m) 

h e Hydraulic head (or englacial water-table elevation) 

(m) 

h s Ice surface elevation (m) 

j Given Julian date (day) 

j c Julian date of melt cessation (day) 

j Q Julian date of melt onset (day) 

tin Conduit meltwater production (kgm -1 s -1 ) 

n Glen law exponent 

n c Conduits per meter in the cross-flow direction 
(rrT 1 ) 

r Conduit radius (m) 

fmax Maximum conduit radius (m) 

r min Minimum conduit radius (m) 

u b Basal sliding velocity (ma -1 ) 

u s Ice surface velocity (ma -1 ) 

u Depth-averaged ice velocity (ma -1 ) 

w Cross-flow width (m) 

x Distance upstream from terminus (m) 
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